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Abstract 

We review the status of quantising (non-abelian) gauge theories using different 
versions of a Hamiltonian formulation corresponding to Dirac's instant and front 
form of dynamics, respectively. In order to control infrared divergences we work in 
a finite spatial volume, chosing a torus geometry for convenience. We focus on the 
determination of the physical configuration space of gauge invariant variables via 
gauge fixing. This naturally leads us to the issue of the Gribov problem. We discuss 
it for different gauge choices, in particular finite volume modifications of the axial 
gauge. Conventional and light-front quantisation are compared and the differences 
pointed out. 
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1 Introduction 



The quantisation of non-abelian gauge fields |]J is notoriously difficult. The reason for 
this is the gauge symmetry: at the level of, say, the Lagrangian a gauge theory is formu- 
lated in terms of redundant degrees of freedom which change under gauge rotations. The 
Lagrangian itself is invariant such that the original and the rotated configurations corre- 
spond to the same physics. Thus they should be identified. In a slightly more formal way 
this can be formulated as follows: Denote by A the total configuration space of the gauge 
theory, i.e. the space of all gauge fields A G A. The space A is unphysical, because it is 
"too big" : it contains an infinite number of gauge equivalent, hence physically identical, 
configurations. If U(x) G Q is an element of the gauge group Q, and the gauge potential 
A & A, then a gauge equivalent field is obtained by the gauge transformation 

A — >• U~ l AU - ~U~ x dU = U A , (1.1) 
9 

where g denotes the coupling constant. In order to get rid of the gauge copies of the field 
A, one identifies all of them to form a single object (equivalence class) called the orbit 
O a of A, 

O a = { u A:U EG} . (1.2) 

A gauge transformation thus leads from one point on an orbit to another point on the 
same orbit, so, by construction, it can never leave an orbit. Two gauge fields on the same 
orbit are called gauge equivalent. The set of all orbits is the physical configuration space 
A4. In mathematical terms this can be written as 



M = A/G, (1.3) 

thus M. is the space where all gauge equivalent fields have been "divided out". Loosely 
speaking one has the identity^ 

A = MxQ, (1.4) 

which states that the total configuration space A consists of the gauge invariant, physical 
variables forming M. = A/g and the gauge variant or "cyclic" variables representing the 
gauge group g. The main problem with gauge theories is to disentangle these two types of 

1 Formula (|l.4[) does not hold in a mathematically correct sense because it says that A is a trivial 
^-bundle over A4. That this is not generally true is the mathematical content of the Gribov problem (see 
below) . 
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variables and make the gauge symmetry explicit by eliminating the redundant variables. 
The resulting theory would then be formulated in terms of the gauge invariant degrees of 
freedom only. This is very much like formulating a rotationally invariant system in terms 
of the radial variable only, after eliminating the cyclic dependence on the angles. In a 
gauge theory, the problem is to properly divide the variables into "radial" and "angular" 
ones. 

Technically, the existence of redundant degrees of freedom shows up through the ap- 
pearance of constraints, i.e. non-dynamical equations between the field variables J|. In 
the Dirac formulation for constrained systems 0, ||, one finds at the first stage the van- 
ishing of the momentum conjugate to the field component Aq (the "primary constraint"). 
However, the real, physical constraint is not this but a secondary one, Gauss's law, which 
is intimately connected with gauge invariance as will be clarified later. It is this constraint 
that has to be solved in one or the other way. In other words, it has to be guaranteed 
that Gauss's law holds once and forever. Then gauge invariance is essentially preserved. 

The formal expression (|1.3|) is of no practical use for finding the physical configuration 
space M.. There are basically two methods to deal with the abundance of variables in 
a gauge theory and "divide out" the gauge group. In the first, one tries to construct 
explicitly gauge invariant variables out of the original gauge fields and formulate the 
whole gauge theory in terms of only these. This was pioneered by Mandelstam both for 
QED |5|] and QCD ||. In this approach one trades the gauge potentials for Wilson loop 
variables |F| which are traced, path ordered exponentials, 

W 7 [A] = -^trPexp j> dx^A^x) , (1.5) 

with 7 a closed path in space-time, P the path-ordering operator along 7 and N the 
dimension of the representation of the gauge group. The gauge theory is thus formulated 
in "loop space" , which, unfortunately, still contains redundant degrees of freedom, as can 
be seen from the appearance of the "Mandelstam constraints" |], ||. Only recently there 
has been some progress in solving these under special circumstances 0. 

The second approach, which will be persued here, aims at identifying A4 via gauge 
fixing, i.e. by finding a subset of A that is isomorphic to A4 as explained e.g. in [10| . 



Such a subset is called a fundamental modular region (FMR) |lTJ- In order to find it, 
one imposes conditions on the gauge potentials A such that one selects a representative 
A on each orbit. Ideally, one would like to find a gauge condition C[A] = 0, such that the 
physical configuration space M. can be identified with the subset 

r = {A e A : C[A] = 0} . (1.6) 
Two requirements have to be fulfilled by an admissible gauge condition : 
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(I) Existence: 

It should select a representative on any orbit; this means that for any A G A there should 
be a solution U G Q of the equation 



A = U A , (1.7) 

with the representative A obeying the gauge condition C[A] = 0. Then any gauge field A 
can be transformed to this gauge. 
(II) Uniqueness: 

There should be only one representative obeying the gauge condition on each orbit. If, 
on the other hand, there are (at least) two gauge equivalent fields, Ai, A 2 , satisfying the 
gauge condition, the gauge is said to be not completely fixed. There is a "residual gauge 
freedom" given by the gauge transformation V connecting the copies, 



A 2 = V A X . (1.8) 
It turns out that in the non-abelian case conditions (I) and (II) cannot be satisfied using 



a simple equation C[A] = for gauge fixing fT2|, Thus, in order to eliminate the copies 
and uniquely specify the gauge, one has to choose appropriate conditions in addition to 
the original one. Generically, the former cannot be written in terms of equations but are 
of non-holonomic nature (like e.g. inequalities). It is obvious that the fields satisfying 
(I) and (II) form a subset A C T which constitutes the FMD to be identified with the 
physical configuration space M.. So we have the relations 



A^ = Acrci. (1.9) 

In the abelian case, i.e. for QED, the gauge fixing program is straightforward and one 
can find M. in terms of (|1.6| ). The Coulomb gauge, V • A = 0, for example, eliminates the 
longitudinal photons and one is left with the physical, transverse ones. These correspond 
to the gauge invariant, "radial" variables. As stated above, for non-abelian gauge fields, 
the situation is much more involved. It was first noted by Gribov in 1978, that the 
Coulomb condition does not fix the gauge uniquely [I2|, a fact that has been named the 
"Gribov problem" thereafter. Gribov's result has been generalised to a wide number of 
gauges by Singer |TB[ . Although there has been a lot of work done on this subject since 
then, (see [1 1] and references therein), it is fair to say that the problem has still not been 
quite resolved. This is particularly true for gauge choices other than the Coulomb gauge. 

So far we have not distinguished between the classical and the quantum theory as the 
remarks above apply to both. The quantisation of non-abelian gauge theories is a problem 
in itself, the features discussed above, however, play a dominant role. 
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There are several methods of designing a quantum theory for non-abelian gauge fields. 
Yang and Mills in their original work M used a Hamiltonian formulation derived from a 
covariantly gauge fixed Lagrangian a la QED, thus quantising all four components of the 
gauge potential A^. The authors, however, did not address the question of finding the 
physical states, which in the covariant (or Lorentz) gauge is non-trivial as known from the 
Gupta-Bleuler formulation of QED ||14|| . This has been criticised later by Schwinger who 
suggested to use the Coulomb gauge instead fl5|] . By insisting on relativistic invariance, 
he obtained a Hamiltonian with a rather non-standard form for the kinetic energy. The 
whole expression looked rather awkward for practical purposes. Thus, in 1976, Feynman 
stated that, "Schwinger, after a lot of hard work found the Hamiltonian..., but... progress 



in this direction ceased" |T£|. In the meantime, people had abandoned the Hamiltonian 
formulation in favour of an alternative method, the path integral. The breakthrough 
had been achieved with the work of Faddeev and Popov ]17||, who inserted their famous 
determinant in the path integral measure leading to the appearance of ghost particles. 
These served to guarantee the unitarity of Feynman amplitudes, explaining earlier work 



of Feynman |TJ| and De Witt ||19fl . After that, path integrals were considered to be the 
method of quantising non-abelian gauge fields, and the shortest route to set up Feynman 
rules and perturbation theory. 

With the work of Gribov ||12|| , however, it became clear that the Faddeev- Popov method 
must be incomplete beyond perturbation theory. This lead to a partial renaissance of the 
Hamiltonian formulation, which clearly is not very well suited for perturbation theory, 
but may offer a better understanding of (some) non-perturbative phenomena. This article 
comments on the present state of affairs and discusses some new directions. 

It is organised as follows: In Section 2, we formulate gauge theories (including the 
Gribov problem) in the simple setting of quantum mechanics. Non-abelian gauge fields 
are introduced in Section 3, where we shortly touch upon the Coulomb gauge and then 
proceed with axial gauges in more detail. In this section we use the conventional equal- 
time quantisation based on Dirac's instant form of relativistic dynamics. The last section 
is devoted to what is called light-cone quantisation (originating from Dirac's front form) 
applied to Yang-Mills theory. 
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2 Quantum Mechanical Gauge Systems 



The gauge theories usually considered are (local) field theories, which genuinely have an 
infinite number of degrees of freedom!. The aspect of gauge symmetry, however, is not 
a field theoretic feature and can be equally well discussed for systems with only a few 
degrees of freedom. Many of the general problems can be discussed there in a simpler 



setting [22, E3, B 



Consider for example a free (point) particle moving on a plane [ 22], |23 1 . Its configuration 
space is 

S={(?i,« 2 )6R 2 }, (2.1) 

and its Lagrangian 

2 2 

L( P ,q)=Piqi+P2q2-Y- J- ( 2 - 2 ) 
From the kinetic term we can read off the canonical Poisson brackets, 



{q ij p j } = 6 ij , i,j = 1,2 (2.3) 

As an additional condition we impose that the angular momentum of the particle should 
be vanishing, 

G = qxp 2 - q 2 pi = , (2.4) 

a constraint which we call "Gauss's law" for reasons to become clear in a moment. The 
angular momentum G is the generator of rotations around the origin of the plane, and we 
interpret it as the generator of time-independent gauge transformations so that the gauge 



group is 0(2). Any motion of the particle violating condition (|2.4 ) is thus gauge- variant, 
i.e. unphysical. Equations ( |2.2|) and Q2.4|) completely define our mechanical gauge system. 

On the classical level, the gauge invariance of the Hamiltonian H = p 2 /2 is expressed 
by the vanishing of the Poisson bracket 



{G,H} = 0, (2.5) 

which simultaneously states the invariance of the Hamiltonian and the time independence 
of the angular momentum. In this very simple example it is clear that the orbits are 
circles around the origin of radius 

2 An exception from this rule are topological field theories |2(J, among them Yang-Mills theory in 
d = 1 + 1 defined on a compact manifold p3] . 
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r = \fql + q 2 2 , (2.6) 

the points of which are transformed into each other via gauge rotations by some angle 0. 
Note that the origin is an exceptional point. There, the orbits degenerate into a single 
point which is not transformed under the action of the gauge group. The origin is thus a 
fixed point and the action of the group is not free. The orbit corresponding to the origin 
has dimension zero whereas all others have dimension one. 

In this example, the gauge invariant variable is obviously given by the radius r, as the 
radial motion is associated with vanishing angular momentum, whereas the angle 

= tan- 1 (q 2 /q 1 ) (2.7) 

is gauge- variant or "cyclic". This can be explicitly checked by calculating the Poisson 
brackets 



{r,G} = 0, {0,G} = 1. (2.8) 

Thus, the angle and the angular momentum G are canonically conjugate. As r is the 
gauge invariant variable, the physical configuration space is 

M = K 2 /0(2) = M+ = {r E R : r > 0} . (2.9) 

The quantisation is straightforward. The Hilbert space of states consists of wave functions 
^(r) being independent of which expresses their gauge invariance. The Hamiltonian H 
acting upon these is just the radial part of the two-dimensional Laplacian, 

H = ---^rl. (2.10) 
2 r Or or 

Let us now try to obtain the same results via gauge fixing. The reader should imagine a 
situation where it is too hard to find the gauge invariant variables, or the ones found are 
too awkward to use. In this case, gauge fixing offers a pedestrian's way to the physical 
configuration space. 

For the case at hand, we choose the gauge 

x(<Zi,?2) = <Z2 = 0. (2.11) 

This gauge exists, since for any point with coordinates (qi , q 2 ) one finds a gauge rotation 
U[4>] transforming it to a point (?i,0) with vanishing q 2 , i.e. on the horizontal axis. The 
angle of rotation is given by = tan _1 (g , 2 /g , i). However, the gauge fixing is not unique, 
since for any point (<?i,0) on the gi-axis there is a copy (— q~i,0) obtained via rotation 
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by 7r or, equivalently, a reflection at the origin. Accordingly, there is a residual gauge 
invariance, not fixed by the gauge choice ( |2.11| ), given by the reflections q± — > —q\. 

This observation can be cast in a more mathematical frame. Consider the Faddeev- 
Popov (FP) "matrix" [T7| 

FP = { X ,G} = q 1 = FP(q 1 ), (2.12) 

which in this simple example degenerates to a single real number. Thus, its modulus is 
the FP determinant, 



J= |detFP| = . (2.13) 

Both quantities, however, are not constant, but dependent on the coordinate q\. The point 
where the determinant is vanishing, qi — 0, is seen to be a Gribov horizon separating the 
two gauge equivalent regions qi > and q\ < 0. For a complete gauge fixing one has to 
demand in addition that qi should be nonnegative, qi > 0. This non-holonomic condition 
defines the FMD of the problem. Again we have consistently obtained M. = K + , and we 



can identify q\ with the radius r, which follows from the definition (|2.6|) upon inserting 



the gauge condition (|2.11|) . This leads to another interpretation of the FP determinant: 



it is the Jacobian of the transformation from the original, cartesian coordinates (qi, ^2) to 



the curvilinear polar coordinates (r, 0) |22 , where, after complete gauge fixing, the radial 
variable r can be identified with the cartesian variable q\ restricted to the FMD. 

Let us discuss the gauge fixing procedure in the quantum theory in more detail. Before 
gauge fixing, the quantum Hamiltonian is 

H= 1 -(pI+P 2 2 ), (2.14) 

with the canonical commutators 

[quPj] = i5ij . (2.15) 

In order not to get a contradiction to the commutation relations we can only impose 
Gauss's law weakly, i.e. on the (physical) states, 

G\^) = 0. (2.16) 
The Schrodinger wave functions are defined as usual, 

*(?!,<&) = <?!,<&!*>, (2-17) 
The momentum operators acting upon them are 
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This is the usual Schrodinger quantisation for cartesian coordinates. Gauge fixing in the 
quantum theory is done by demanding that the physical states (satisfying Gauss's law) 
are given by a projection onto q 2 = 0, 

* x (gx)= (ft,fc|*>| B=a0 = <fc,0|¥>. (2.19) 

On the gauge fixed physical states the momentum p 1 acts as the usual differential operator, 
Pi = —id/dqi, but, due to Gauss's law, 

Pal*) = ftWl*) 7^ -*^-|*) • (2-20) 



dq 



2 



This yields pi 



(gi,0|pl|*) = (gi,0|p 2 (?rW)l*) 



and thus the Hamiltonian 



--^-*xfe) , (2-21) 
<?i oqi 



acting on gauge invariant states. As expected, it coincides with (|2.10| ) upon identifying 
J(Qi) = with r. The Hamiltonian ( |2.22|) is hermitean with respect to the scalar 
product 

(*i$> = ^ j d 2 q^{q)^{q) = ^ dq l q 1 %(q 1 )^ x (q l ) , (2.23) 

where we have divided by the "volume of the gauge group", 2n. The associated Hilbert 
space of physical states is thus 

H = L 2 (R+,J(q 1 )dq 1 ) . (2.24) 

The global properties of the configuration space M.q, which is a half-space with a boundary, 
affect the wave functions in a peculiar way. One can imagine the configuration space as 
being a semi-infinite potential well, so that the wave function can not penetrate the barrier 
to negative values of q\ or r. This necessitates the boundary condition 

* x (qi = 0) = . (2.25) 
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The wave function is thus repelled from the origin. Another way to see this is to simplify 
the measure in the scalar product by absorbing a factor J 1 / 2 into the wavefunction. This 
results in a rescaled Hamiltonian 

H = J-^HJ 1 ' 2 = -~ - , (2.26) 

which has a standard form for the kinetic energy but an additional "effective potential" 
term l/8q 2 , reminiscent of the centrifugal barrier in the hydrogen atom. The new Hamil- 
tonian is hermitean on L 2 (Rj , dq\ ) . 

Let us close this section with a mathematical remark. It should be noted that one 
cannot write the configuration space M 2 as a trivial bundle over i.e. as a product of 
the physical configuration space with the gauge group, 

M 2 ^ 0(2) x M+ , (2.27) 

because of the singularity at the fixed point r = 0. This is tantamount to saying that the 
mapping from cartesian to polar coordinates, 

{qi,q 2 )^{r,<t>) , (2.28) 

is not one-to-one. At the origin, r = 0, the associated angle <fi is not defined. This 
phenomenon, the Gribov problem, does equally occur within Yang-Mills theory, which is 
our next topic. 



3 Hamiltonian Yang-Mills Theory: Instant Form 



3.1 General Remarks 



In the instant form of relativistic dynamics |25| one choses the ordinary Cartesian time t = 
x° as the parameter describing the causal evolution of a given system. The Hamiltonian 
then propagates a state according to the Schrodinger equation, 

HW) = i^Wi- (3-1) 

This choice of time is most common and coincides with the non-relativistic one, which is 
unique due to Galileian causality. Einstein causality, however, leaves room for alternative 
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choices, one of which, leading to the front-form of dynamics, will be discussed in the next 
section. 

First of all, let us explain our basic conventions and notations for gauge field theory. 
Throughout this paper we will use a finite spatial volume as an infrared regulator, following 
the approach pioneered by 't Hooft ]2E| and Liischer P7| . We let — L < < L, i = 1, 2, 3, 
and impose periodic boundary conditions (pBC) on the gauge fields, 



-L) = Ai{xj = L) , i,j = 1,2,3. (3.2) 
This results in a torus geometry for three-space. We write the YM Lagrangian in the 



following form 28, 29 



= EtAZ-±(EtEf + BtBt)+AtD?E$, (3.3) 

where Af denotes the time derivative of the field (the velocity), H the Hamiltonian and 
G a the Gauss operator. For simplicity, we choose SU{2) as the gauge group and will 
switch freely between component and matrix notation for the gauge fields, Aj = A°;T a , 
the chromo-electric fields, 

Ei = E«T a = A t - Aio , (3.4) 

and the chromo- magnetic fields 

Bi = B?T a = e ijk (dj + igA^Ak . (3.5) 

As a basis for the Lie algebra su(2) we have chosen T a = r a /2, with r a the usual Pauli 
matrices. The covariant derivative is 

Df = + ge abc A1 , (3.6) 

or in su(2) matrix notation, 

Di = di + ig[Ai, } = d { + ig adA • (3.7) 



From the kinetic term of (R.3I) we read off 29 the canonical commutators 



[^(x,t),^(y,t)] = dy«(x-y) . (3.8) 

Thus, gauge fields A { and chromo-electric fields Ei are canonically conjugate and generalise 
the cartesian coordinates q, p of Section 2. It is therefore straightforward to set up a 
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functional Schrodinger picture with wave functionals ^[Ai] replacing the wave functions 
of quantum mechanics. The canonical momenta act as 



EMM = -^nM , (3.9) 
from which one can infer the commutation properties of the Gauss operator 

[G a (x), G\y)} = ^ c G c (x)5(x - y) , (3.10) 

which yields a representation of the Lie algebra su(2). As emphasized by Jackiw ||30|| , the 
field theory discussed so far might be totally consistent but differs from the gauge field 
theory defined by additionally imposing Gauss's law on the (physical) states, 

G^IA,} = . (3.11) 

The Gauss operator generates time-independent, infinitesimal gauge transformations ac- 
cording to 

[A°, G b ] = iDf = iS ab di + ige ahc A\ , (3.12) 

where we have omitted the spatial arguments for convenience. The first term on the r.h.s. 
of ( |3.12|) corresponds to a translation in ordinary space, the second, non-abelian one, 
to a rotation in color space. Chromo-electric and -magnetic fields, however, transform 
homogeneously, 



[E^G b ] = ige abc E^ (3.13) 

[Bf, G b ] = ige abc B^ . (3.14) 

The Hamiltonian is (gauge-) invariant under ( |3.12| ), 

[H,G a }=0 } (3.15) 

thus Gauss's law is consistent with the equations of motion. Finite gauge transformations 
are generated by the operator 

Sl[6] = exp -i J d 3 xE?(x)Df6 b (x) , (3.16) 
such that the transformed field is given by 

nAitt~ l = U- l A t U - -U~%U = % , (3.17) 

9 

where U = U[9] is defined as 
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C/(x) = expig9 a (x)T a G 5C/(2) 
The finite version of Gauss's law, ( 3.11] ), is then 



(3.18) 



n%) = *[AJ , (3.19) 

expressing the gauge invariance of the physical states. Let us now try to find the physical 
configuration space via gauge fixing. Before we discuss any particular example, some 
general remarks are in order. We denote the gauge fixing condition by 

X [A t ] = (3.20) 
and imagine a coordinate transformation 



A, 

n, 



(A*,x) 

(n*,7r) 



(3.21) 
(3.22) 



where IL = Ei is canonically conjugate to Ai, and n to \. The gauge invariant canonical 
pair is (A*, II*). The ideal case seems to be that tc coincides with the Gauss operator G, 
and that the gauge condition \ is conjugate to it. Gauss's law would then state that the 
wave functionals are independent of x (being the analogue of the angle of Section 2), 



G^[Al x ] = - i —J[Al x ] = 



(3.23) 



and we would have found the physical configuration space Ai consisting of the A*. In this 
case the FP matrix is (in symbolic notation) FP = G} = 1, its determinant has no 
zeros and the gauge fixing were complete. Following Gribov's original observations 
Singer has proven that such a scenario cannot possibly be achieved (T3| within gauges 
allowing for one-point compactification, IR 4 U {oo} = S 4 |HJ, so that the situation is 
considerably more complicated. This can be easily seen [23| from (|3.23|) : if the Gauss 
operator could really be represented as the (functional) differential operator of (|3.23|) , it 
would commute with itself, in contradiction to ( |3.10|) . Clearly, no such contradiction is 



present in the abelian case. Similar arguments have been given recently in [32 



Let us consider the FP matrix for arbitrary gauge fixing \- m terms of the original 
variables A one finds 



FP = {x,G} 



5Ai 



D, 



(3.24) 



So, if x is linear in the Ai, the FP matrix is essentially the covariant derivative. In terms 



of the new, gauge fixed variables A*, %, the FP matrix is [28 
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FP = ^ . (3.25) 

OTT 

From this last expression one infers that if the FP determinant is non-vanishing, 

J = | det FP| = | det — | 7^ , (3.26) 
on 

Gauss's law can be solved (locally) for the gauge variant momenta it by inverting the 
FP matrix. If this can be done explicitly, one is left with the gauge invariant variables 
(A*,U*). It will turn out, however, that also this program is beset with obstructions. In 
general, the FP matrix is a (partial) differential operator. Ideally, one would like to know 
its whole spectrum. Anyhow, in order to invert it, one has to specify its domain together 
with boundary conditions. This is another, if technical, reason to work in a finite, definite 
geometry like the torus. Still, the FP matrix might be too complicated an operator to be 
inverted exactly. 

It should again be noted that, although the FP matrix was defined within the classical 
theory in (|3.24|) and ( |3.25|) , it has a strong impact on the quantum theory. Using a 



functional Schrodinger picture, a FP matrix that is only dependent on the Aj, (which is 
the generic case), acts as a pure number on the functionals ^LAJ. The remarks above 
thus equally apply to the quantum theory, when Gauss's law is solved in the course of 
defining the quantum Hamiltonian acting on the physical states via 

H X * X [A*] = |*> x=0 . (3.27) 
Let us now discuss some prominent examples of gauge fixing. 



3.2 The Coulomb Gauge 



The Coulomb gauge condition, which is the most natural one in the abelian case 
demands transversality of the gauge potentials, 

x[Ai) = diAi = . (3.28) 

Within the Hamiltonian formulation, it is the most widely used and therefore best stud- 
ied gauge fixing condition (see for example |Tl|, |l^, ^2[ [34], |35]). The existence has been 



implicitly proved by Semenov-Tyan-Shanskii and Franke [Q. The gauge invariant vari- 
ables are (like in QED) chosen to be the transverse gauge potentials, A* = Af , such 



14 



that the physical states are \1/ = The canonical momentum, the electric field, is 

decomposed into a longitudinal and transverse part JT5]], 

Ei = Ej + Ef = Ej + d4 , (3.29) 

where a scalar potential </> has been introduced. The latter is canonically conjugate to the 
gauge condition, 

{ X ,0} = 1, (3.30) 

therefore corresponds to the gauge variant momentum 7r and should be eliminated via 
Gauss's law, (to be imposed on the physical states), 

G\ x=0 = diDi[A^]<f> + [Aj, Ef] = FP<fi — p c = . (3.31) 

This expression might be called the "non-abelian Poisson equation" with the Laplacian 
replaced by the FP operator, FP = diD i} and the color-charge density of the gluons, 
p c . As stated above, it can be solved for </> if the FP determinant J = | detFP| is non- 
vanishing. The latter is difficult to check, since the FP matrix is a complicated partial 
differential operator. Its inverse is not explicitly known and can only be approximated as 
a power series in the coupling g. Formally this means 



FP" 



0(g), (3.32) 

so that the perturbative expansion can be viewed as an expansion around the abelian 
limit. As a result, the Hamiltonian, which depends on FP -1 and J, can only be written 
down in a formal or a perturbative sense [|15[ | 



As mentioned above, Gribov was the first to point out |12[ that for large fields, i. e. be- 
yond perturbation theory, the FP matrix develops vanishing eigenvalues on what are now 
called "Gribov horizons". The (convex) interior of the first horizon, the Gribov region B, 
however, still contains Gribov copies and has to be further restricted to form a FMD. This 
is most transparently discussed in terms of a (Morse) functional on A. Its critical points 
are the fields A{ satisfying the Coulomb gauge condition, its relative minima constitute 



the Gribov region B : its absolute minima the FMD C B [10, [11], |34|| . 

In view of the above, in particular ( |3.32|) , one can say that the Coulomb gauge is well 
suited for perturbation theory ||32|| . This has been used e.g. in the Hamiltonian proofs 



of asymptotic freedom |36| or in Liischer's small volume calculations |p7 |. To go beyond 



perturbation theory is possible but technically involved |TT , |37j . It is therefore desireable 
to look for and find alternatives to the Coulomb gauge formulation. 
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3.3 Axial- Type Gauges 



Axial gauges have been introduced independently by Kummer |35|] and Arnowitt and 
Fickler |3P| in the early sixties (for recent reviews see |4T], They are defined by the 
condition 



n-A = 



(3.33) 



where n is a constant four vector. Depending on whether n is time-like, space-like or 
light-like, one refers to the temporal, the axial or the light-cone gauge, respectively. Their 
algebraic simplicity was - and still is - the main motivation for using them, despite their 
lack of explicit covariance. As the FP operator, FP = n ■ d, is field- independent, the 
gauge is (apparently) ghost free pE2| . Accordingly, the constraint equations are rather 
easy to solve, one does not have to transform to curvilinear coordinates and there are no 
ordering problems, in contrast to the Coulomb gauge. But it has already been noted by 



Schwinger in 1963 |43| that there is a huge residual gauge freedom: gauge transformations 
independent of n • x are compatible with the gauge condition. These transformations 
correspond to zero modes (ZMs) of the FP operator, thus, det FP = det(n ■ d) = 0. 
As pointed out by Schwinger, this leads to infrared divergences in the Hamiltonian. The 
solution to the infrared problems is again to work in a finite spatial volume, — L < Xj < L, 
and impose pBC (torus geometry). There is however a price for this to be paid: the (strict) 
axial gauge does not exist on a torus. For the abelian case, this has been first noted by 



Yao []44], a student of Schwinger's. For the non-abelian case, it is implicit in the work of 
Batakis and Lazarides [^5| on the vacuum of finite-temperature Yang-Mills theory, where 
pBC naturally arise in the (imaginary) time direction so that the temporal gauge cannot 



be attained (see also ^6|). For the light-cone gauge, the same was first noted by Franke 
et al. j47|, for space-like axial gauges it is implicit in fH) and explicitly stated in 

mug, HI- 

Let us for the sake of definiteness consider the (space-like) axial gauge, A 3 = 0. It 
turns out then, that configurations a 3 (xi,x 2 ), which are independent of x 3 and thus ZMs 
of FP = <9 3 , cannot be gauged away by allowed gauge transformations!. Another way to 
see this is the following. Consider the Wilson loop winding around the torus, 



1 f 
W[A 3 ] = —trP 3 expig J dx 3 A 3 , 



(3.34) 



3 These are periodic up to an element of the center of the gauge group and thus preserve periodicity 
of the gauge fields p7|. The center of SU(2), for example, is Z 2 = {±1}- 
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where P 3 denotes path ordering in x 3 direction. The Wilson loop is a gauge- invariant, 
dynamical quantity. In general, therefore, ^[^3] 7^ 1, whereas W- 7 ^ = 0] = 1. If we 
decompose A 3 into a ZM 03 and a non-ZM A 33 , 

1 f 

A 3 = A 33 + a 3 , a 3 = — J dx 3 A 3 , (3.35) 

— L 

one can only gauge away A 33 but has to retain a ZM a' 3 , the gauge transform of 03. This 
corresponds to the modified gauge condition 



d 3 A 3 = . 



(3.36) 



Thus, the Fourier transform A 3 (k) is zero everywhere, except for k 3 = [44]. One should 
also be careful to realize that in the non-abelian case, 



1 f 1 

W[A 3 ] = —trP 3 expig J dx 3 A 3 ^ — tr exp[2^La 3 ] 



(3.37) 



due to the path-ordering operation. Also, with the gauge condition (|3.36|) , one still has 
to fix the residual gauge freedom of gauge transformations independent of x 3 . There are 
many possible choices [44|, [49], |50|, |52|, [53|. We will follow Palumbo |0| [51]], and call this 
choice henceforth Palumbo gauge. To define it, we introduce a momentum lattice via 
Fourier decomposition into plane waves, 



e(n ■ x) = exp in(nix 1 + n 2 x 2 + n 3 x 3 )/L 
For an arbitrary phase space variable, F, we write 



(3.38) 



F(xi,x 2 ,x 3 ) = F(n 1 ,n 2 ,n 3 )e(nix 1 +n 2 x 2 + n 3 x 3 ) 

n.1,712,713 

= F(n 1 ,n 2 ,n 3 )e(n 1 x 1 + n 2 x 2 + n 3 x 3 ) + 

n 3 #0 

+ Y F( n i, n 2,n 3 = 0)e(niXi + n 2 x 2 ) + 



n 2 #0 



+ Y F ( n ii n 2 = 0, n 3 = 0)e(n 1 x 1 ) + 

+ F(m = 0, n 2 = 0, n 3 = 0) 

= F 3 (xi,x 2 ,x 3 ) + F 2 (xi,x 2 ) + Fi(xi) + F . 



(3.39) 



In this way we have split the function F into four components F r depending on less and 
and less space coordinates Xi, 



17 



F(x) = E F r . (3.40) 

r=0 



In coordinate space, the components F r read explicitly 



1 ) 

F 3 (x 1 ,x 2 ,x 3 ) = F(x) - — J dx 3 F{x) , (3.41) 



—L 

L 



F 2 (x 1 ,x 2 ) = J dx 3 F(x) - — J dx 3 dx 2 F(x) , (3.42) 



(2L)= 



1 L 1 L 

Fl(;ri) = J^l L dX ' dX2F{X) ~J^lf XF{X) > (3 ' 43) 



l r 

F ° = mysj d * xF ^ ■ ( 3 - 44 ) 



All phase space variables (gauge fields and momenta) will be decomposed analogously. 
In the following, all summations over r will be written explicitly to avoid confusion. The 
components A ir , r > 0, do not contain a ZM with respect to x r , 

L 

J dx r A ir = , r > , (3.45) 

— L 

and satisfy 

d r A is = , r > s . (3.46) 
The Palumbo gauge is then defined by the following conditions 



^33(^1,^2,3:3) = 0, 
A 22 (xi,x 2 ) = 0, 

A 11 (x 1 ) = 0, (3.47) 

which in shorthand notation is simply 

A rr = 0, r = l,2,3. (3.48) 



A last definition we need is 



r-l 

Gr = E , (3.49) 

s=0 
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representing the sum of all ZMs with respect to x r contained in A r , thus 



drdr = . 



(3.50) 



Our next task is to show that the Palumbo gauge (|3.48|) exists. This proof has not been 



given in the original work of Palumbo [49], |51|, who simply assumed the existence. As a 
consistency check he solved Gauss's law for the momenta conjugate to the A rr without 
redundant degrees of freedom being left. One can, however, give a direct existence proof 
]55| by constructing a gauge transformation U that takes an arbitrary configuration A to 
Palumbo gauge, with 



u A rr = OJ- l A r U - -U^drU), 
9 







(3.51) 



To do so, one has to perform three steps for r = 1, 2, 3, such that U = UiU 2 U 3 , d r U s = 
for r > s. The following table schematically depicts the procedure. 



before gauge fixing 


after gauge fixing 


A 3 = A 33 + a 3 


A' 3 = + 4 


A 2 = A 23 + A 22 + a 2 


A' 2 = A' 23 + + a' 2 


A 1 = A 13 + A n + A n + a x 


A[ = A' 13 + A' 12 + + a[ 



Table 1: The gauge transformations to Palumbo gauge. In the right column, the prime 
denotes successive application of the transformations U 3 , U 2 and U\, respectively. 

The explicit form of the gauge transformations U r , r = 1,2, 3, is 



P r exp 



-w 



exp [ig(x r + L)a! r 



= h r [x r ; A r ] exp [ig(x r + L)a' r ] 



(3.52) 



The first, path-ordered exponential, h r , gauges away "too much", namely A rr and the ZM 
a r and therefore cannot be (and is not) periodic in x r . Periodicity and a ZM a' r are then 
restored by the second exponential [[47], |52| . It should be noted that the ZM a' r differs 
from the original one, a r due to (|3.37| ). One has 



2gL 



lnh[x r 



L; A r ] , 



(3.53) 



which is, however, formal in so far as the logarithm on the r.h.s. is a multivalued func- 
tion. This raises the question of uniqueness which will be adressed shortly when we have 
calculated the FP determinant. 



19 



Before we come to that, let us write down Gauss's law in our decomposed notation. The 
spatially constant part (r=0) will be discussed elsewhere (see also Palumbo's discussion 
P§| , |5"T|]). For r > 0, Gauss's law reads in components 



G r = d r E rr = r E rr + ig[a r , E rr \ = FP^a^-E^ = p r . (3.54) 

This should be solved for the components E rr , which are canonically conjugate (in the 
Poisson bracket sense) to the gauge conditions \r = A„. To this end, one has to invert 
the FP operator j54j, 



FP[a r ] 



(3.55) 



/ d^ax] \ 

ig &dA 12 d 2 [a 2 ] 
^ ig &dA 13 ig adA 23 d 3 [a 3 ] ) 

where the operation 'ad' has been defined in ( |3.7| ). The FP operator is field dependent and 
in a path integral formulation a la Faddeev and Popov would lead to ghosts. Nevertheless, 
it is still much simpler than the Coulomb gauge expression, because it is an ordinary rather 
than a partial differential operator (for each r). Due to the algebraic property, 



d r a r = = [d r , i 



(3.56) 



the eigenvalue problem for d r furthermore separates into space and color. As a result, 
it is exactly (i.e. non-perturbatively) solvable |50| , pT , |52| . The eigenvalue problem (for 
SU(2)) is given by the expression 



FP ab ui 



(d r 5 ab + ga c r e a c >L ar = K r , ar 



,c ,ab\ . ,b 



The eigenf unctions are found to be 



Mr 



exp(i7in r x r / L)w ar 



and separate into a plane wave times a color matrix w. The eigenvalues are 



A, 



irn r 



n r 



±1,±2,...^0 



a r e 



|o, ±i|a r || 



(3.57) 



(3.56 



(3.59) 



where \a r \ = (a a r a a r ) l/2 is the length of the color vector with components a r . 

In order to dicuss the uniqueness of the gauge fixing, we calculate the FP determinant 



J r = detFP(a r )/detFP(0) = IJ V,a,/ II K,a r =0 = 

= ((7|a r |L)- 2 sin 2 (^|a r |L) . (3.60) 
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It is thus essentially given by the Haar measure of SU(2) pl| , |5"1"| , [5^ ] and vanishes for 
configurations a r of length 



\a r \ — im r /gL 



(3.61) 



which are gauge equivalent to a r = 0. Thus one has an infinity of Gribov copies labelled 
by the integers n r . The fundamental modular region is therefore given by a "ball" B in 
color space, 



B : < \a r \ < n/gL 



(3.62) 



In the region B, Gauss's law can be solved uniquely with the help of the inverse of FP, a 
Green function 



FP;\x r ,y r )= £ 



U n r ,a r ( x r)u nrj a r (y r ) 



„ r#0 iim r /L + ga r 

Ot-p 

Our Hilbert space of states consists of the functionals 



(3.63) 



H 3 V[A 13 , A 23 , A l2 , a 3 , a 2 , ai] = ^ X [A*} , (3.64) 

where Ai%, A 23 can be viewed as transverse gluons (with the ZMs removed) and all the 
other fields are ZMs or lower- dimensional fields |52|]. The state functionals have to vanish 
at the boundary of the fundamental region B, 



laJ = 



|a r | = 7r/ 1 gL 







(3.65) 



The quantum Hamiltonian is implicitly defined on the gauge invariant states via 



H^ X [A*} = (A\H\*) Arr=0 
and explicitly found to be (integrating over the torus T), 



(3.66) 



5 



S A- a A 

r=o u-tiir u -f^ir r=0 



+ 



j4 7' i 



+ £ ^'(FP; 1 * PJMFP; 1 * Pr)\ 
r=0 ) 



(3.67) 



The first term containing the functional derivatives is the kinetic energy, the term in 
the second line is the analogue of the Coulomb terms stemming from the solution of 
Gauss's law, the asterisk (*) denoting convolution. As was the case for the Coulomb 
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gauge, the appearance of the Jacobians J r guarantees the hermiticity of the Hamiltonian 
in the physical Hilbert space. This in turn ensures the unitarity of the S-matrix. We thus 
note the following rule. If the FP matrix is field dependent, the (quantum) Hamiltonian 
contains non-trivial Jacobian factors. In the path integral formulation there appear ghosts. 
Either of these features is necessary to maintain unitarity of the theory. This adds another 
point of view to the debate whether continuum axial gauges, A 3 = 0, are ghost-free, or 
contain ghosts upon infrared regularisation |^8|, In our opinion, it is still an unsettled 
problem whether there is a consistent (unitary) formulation of axial-gauge QCD that is 
definitely ghost-free. 



4 Hamiltonian Yang-Mills Theory: Front-Form 



In the light-cone (LC) formulation of quantum field theory based on Dirac's front-form of 
relativistic dynamics p5| one introduces new variables x ± = 2 _1 / 2 (x° ± x 3 ) and specifies 
canonical commutators on hyperplanes x + = 0, tangent to the light-cone (for recent 
reviews see [58|). This drastically alters the canonical structure of field theories. They 
typically have Lagrangians linear in the LC velocity; some variables become redundant via 
second class constraints and can in principle be eliminated from the theory. A Lagrangian 
linear in the velocity is singular in the sense of Dirac [|J and might be treated with his 
theory for constrained systems. More appropriate, however, is the method of Faddeev 
and Jackiw [p9fl , which is explicitly taylored for this case. It is essentially equivalent to 
Schwinger's quantum action principle |57|, which will be shown elsewhere [^9(1 . 



For our formulation of LC Yang-Mills theory it is again convenient to work a in finite 
volume, —L < x~,Xi < L, i — 1,2, and impose pBC in every spatial direction. The 
Lagrangian is 



c = n a A a _ + n%A%--(n a n a + B a B a )+A a + (D ab n b + Dfn*) 

= U a A a + U^Al -H + A a + G a , (4.1) 

where 7i denotes the Hamiltonian (density), G a the Gauss operator and the dotted vari- 
ables derivatives with respect to the LC time x + . Furthermore, we have defined the 
chromo-magnetic fields 



B = F 12 = d,A 2 - d 2 A, + g[A u A 2 ) (4.2) 



22 



and the momenta 



IT = A a _- D a *A b + , (4.3) 
= d.A* - DfA b _ . (4.4) 

The Lagrangian ( |4.1| ) appears to be perfectly canonical. However, as is typical for the 
LC formulation, the momenta 11° are dependent variables as ( |4.4j ) is a (second class) con- 
straint, which expresses the IT? in terms of other degrees of freedom. The number of 
degrees of freedom per space-time point is thus only 12 versus 18 in the instant formu- 
lation. Upon inserting (|4.4|) into ( |4.1|) , the Lagrangian becomes highly non-canonical. 
Nonetheless, the method of Faddeev and Jackiw |29| can still be applied and gives the 
following elementary brackets (in condensed notation) 



{Ai,Aj} = ^DZ 1 , (4.5) 

{A,T1} = -DiDZ 1 , (4.6) 

{A.,U} = 1, (4.7) 

{IT, IT} = D i DZ 1 D l . (4.8) 

These brackets look highly non-canonical (i.e. non-cartesian, as they differ from zero or 
one) and furthermore are field dependent. The inverse of the covariant derivative _D„L4_] 
is difficult to obtain before gauge fixing, even classically. Finally, the program of first 
quantising (preferrably in cartesian coordinates) and then fixing the gauge, as was done 
in the instant formulation, seems rather hopeless to be pursued. One does not even know 
how to define a (functional) Schrodinger picture. The only possible way to proceed seems 
to gauge fix before quantisation, as was already done by Franke et al. [[47} . 



It is illuminating to note that all the problems above disappear in the pure continuum 
LC gauge, A_ = 0, which was used successfully in many perturbative applications (see 
and references therein). In this case, the redundant momenta are 

II? = d_A1 , (4.9) 



and the elementary brackets are reduced to the canonical LC bracket 



{Ai,Aj} = \5ijdZ 1 ■ (4.10) 
As in the case of the axial gauge, however, the LC gauge does not exist on a torus. This 



was already noted by Franke et al. j|7| who suggested a modification which we will call 
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the FNP gauge. It is convenient to change to the Cartan-Weyl basis for su{2) matrices 
via 

papa = p+p- + p- T+ + F 3 T 3 ? pa = r « ^ _ 

The FNP gauge is then defined as 

At = 

A 3 _ = v , <9_v = , (4.12) 

where the upper indices refer to color, and the lower ones are spatial. Again, as in the 
Palumbo gauge, a ZM (with respect to x~) is retained which in addition is diagonal in 
color space. (An analogous gauge fixing in the equal-time case was used in [52].) 
We do not use the Palumbo gauge within the front-form dynamics as this amounts to 
treating the indices -,i = 1,2 in a symmetric fashion |54|, |62|| . As the fields A_, A iy after 
elimination of the IT, enter the Lagrangian in a manifestly asymmetric way, this seems 
to be somewhat inappropriate. 

Setting G a = for the moment, the Lagrangian becomes 

£ = nv + At(d_ - igv)A- + Ar(d + + igv)A+ + A*d_A\ - F[n± vr, A t ] , (4.13) 
from which one can read off ^l|] the elementary brackets 



{v,tt} = 1, (4.14) 

{AlA)} = §<VT\ (4.15) 

{Af,AJ} = -\5 ij {d_±igv)- 1 . (4.16) 

These lead in the standard manner to the following FP matrix elements, 



FP±0) = d„±igv, (4.17) 

FP 3 = d_ . (4.18) 

Note that only the first two of them are field dependent (on the ZM v). Gauss's law, 

G = FP*n- DAli = , (4.19) 
can again be solved exactly and uniquely if the FP determinant 

7TTI 

J = | det FP| = [] ("T + a 9 v ) ~ siv- 2 (gLv) (4.20) 
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is non-vanishing. This is guaranteed in the fundamental modular region 



< v < n/gL , (4.21) 

which in the FNP gauge is a "line-segment" in color space. The associated Hilbert space 
in the quantum theory can be found as follows. The off-diagonal components Af have a 
Fock expansion with f-dependent coefficients, whereas the coefficients in the expansion of 
A 3 are f-independent [[|7], |63|]. The canonically conjugate variables v, ir = —i5/5v act on 
a Hilbert space of functionals ty[v] depending on the ZM v. This implies that the total 
Hilbert space is generated by a basis consisting of products of Fock states (associated with 
the transverse gluons Ai) with the functionals ty[v]. The field dependent coefficients in the 
expansion of Af thus are no problem as they act multiplicatively (i.e. like pure numbers) 
in the ^-representation. In summa, this supplements the successfull Fock space picture of 
the continuum LC gauge with ZM degrees of freedom which are generally believed to be 



connected with non-trivial topology ||7, |64], | 65 
For the quantum Hamiltonian one finds 



H= - [ d 3 x\ - J^M-^-JM-^- + B a B a + 
2 J ov ov 

T K 

+ 2FP: 1 * (dtlj+fpt 1 * (An,)- + d: 1 (D l uA 3 d: 1 (D l n l ) 3 ) (4.22) 



where for the redundant momenta one has to substitute the constraint (O) in the 
FNP gauge. This Hamiltonian has been obtained by Franke et al. [47||, who, however, 
overlooked the necessity for the Jacobian J to appear in the kinetic term. Unfortunately, 
the Hamiltonian ( 4.22|) is not the end of story. It still contains redundant variables which 
have nothing to do with a (residual) gauge freedom but are second class constrained and 
analogous to a constrained ZM extensively studied in LC 4 theory [66]. The constraint 
is 



SH_ 

5a? 







where af is the ZM of A 3 with respect to x .It cannot be solved in closed form [47 
the moment, the implications of this constraint are unclear. 



(4.23) 
At 
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